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Abstract 

We consider two bosonic atoms interacting with a short-range potential and trapped in a spher- 
ically symmetric harmonic oscillator. The problem is exactly solvable and is relevant for the study 
of ultra-cold atoms. We show that the energy spectrum is universal, irrespective of the shape of 
the interaction potential, provided its range is much smaller than the oscillator length. 
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I. INTRODUCTION 



From the earliest studies of the dynamics of our solar system to the first attempts to 
apply quantum mechanics to the simplest atoms (hydrogen and helium), few-body systems 
have played a pivotal role. In particular, the problem of two particles interacting via a 
central force law is familiar to students of physics. After elimination of the center-of-mass 
motion the two-body problem is reduced to a one-body problem in the relative co-ordinate. 
The well-known two-body gravitational and the Coulomb problems admit exact analytical 
solutions, which can be tested by experiments. 

The few-body problem is still of great interest, especially in view of recent experimental 
advancements in the area of ultra-cold atoms. In these experiments neutral atoms are 
first cooled to nano-Kelvin temperatures and then confined in an optical lattice formed 
by standing wave laser beams, where there are as few as two-and three atoms per lattice 
site. Typically, the lattice sites confining the atoms may be regarded as harmonic oscillator 
potentials with a trapping frequency u of a. few tens of cycles per second. The coupling 
between the harmonic oscillator wells (that is, between the lattice sites) may be controlled 
by the intensity of the laser beams. In the limit of low tunneling, the individual wells 
may be regarded as independent three-dimensional (3D) harmonic oscillators. The strength 
of the pair-wise interaction between the atoms on each site may be fine-tuned by sweeping 
through a static magnetic field using a Feshbach resonance. ^"^i^ At nano-Kelvin temperatures 
quantum effects dominate, and we thus have what amounts to a non-relativistic quantum 
few-body system whose inter-particle interactions can be controlled at will. This achievement 
is truly remarkable, and such ultra-cold matter experiments can serve as ideal laboratories 
to investigate the quantum few-body problem. 

The experimental group of Stoferle et al? has trapped two atoms in a 3D harmonic 
oscillator well, and by tuning near the Feshbach resonance, can make the inter-particle 
interaction strong enough to form diatomic molecules. Because the two atoms (each of 
which behaves as a fermion^) occupy distinct hyperfine states, they may still interact in the 
relative s-state. Subsequently, the energy of the state is deduced by dissociating the molecule 
using an rf-pulse and measuring the binding energy of the molecule. The binding energy 
as a function of the strength of the interaction (which is related to the s-wave scattering 
length) is mapped out by by sweeping a magnetic field using a Feshbach resonance. In these 
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experiments the atoms form only loosely bound pairs so that the interaction range may be 
taken to be very small compared to the size of the two-body system. This system is then 
an experimental realization of the exactly solvable problem of two harmonically trapped 
particles interacting via a zero-range pair-wise potential at zero temperature. 

In this article we present the exact solution to this quantum two-body problem, assuming 
only a knowledge of quantum mechanics at the senior undergraduate or beginning graduate 
student level. Remarkably, this ostensibly simple problem has only been investigated within 
the last 10 years,-!^ and now enjoys a new life due to its relevance to the experiments 
discussed in Ref. |2l. In particular, the theoretical spectrum of this system is in excellent 
agreement with the experimental spectrum we have described.- Our approach does not make 
any reference to the shape of the short-range two-body potential, and requires only that its 
range be much smaller than the oscillator length,- which serves as a natural measure of the 
particle spacing. This approach beautifully illustrates the notion of shape-independence of 
potentials (for low energy scattering and binding) that was first developed in the context of 
nuclear physics.— Shape-independence is embodied in the effective range expansion, which is 



explained and derived in Ref. lUj. The concepts of scattering length and effective range that 
come in this expansion are important in the analysis of the experimental results with ultra- 
cold atoms, and readers are encouraged to study the excellent discussion of these topics in 



Ref. 
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We shall make use of these concepts and show that the energy levels of this system 
are independent of the shape of the interaction potential in the limit of zero range. This 
result is accomplished by noting that, for a harmonic trap, the center of mass motion may 
be isolated, and the interacting relative motion may be written as a one-body problem. The 
irregular solution of the harmonic oscillator wave function, which is normally discarded in 
the noninteracting problem, will be shown to be the relevant solution for the zero-range 
interacting case. This realization paves the way for a gentle introduction to the concept of 
a pseudopotential in the form of a regularized delta function. 

The zero-range pseudopotential, first introduced by Fermi,— and applied to the many- 
body problem by Huang and Yang,— is extensively used in the current literature of ultra-cold 
gases. The same pseudopotential has also been used to study the spectrum of two particles 
interacting by a zero-range potential in a spherical harmonic oscillator well. Rather than 
taking this path, we shall follow the derivation of Jonsell,^ which is in the same spirit as our 
initial intuitive approach. In Sec. Ill we finish with some concluding remarks and suggestions 
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for future research in this area. 



II. THE TWO-BODY PROBLEM 
A. The zero-range pseudopotential 

We first consider the problem of two identical bosons, each of mass M, in the absence 
of a confining potential. The bosons interact with a short-range central potential. Most of 
what follows in this section also applies to fermionic atoms in the spin-singlet state, or to 
nucleons in the s-state with an antisymmetric spin-isospin wave function.— In the relative 
co-ordinate r = |ri — the s-state asymptotic scattering wave function for positive energy 
E is given by 

u{r) = riplr) ~ sin(A;r + S{k)), (r > b). (1) 

We have E = h^k'^/M, and b is the range of the potential. The effective range expansion, 
relating the phase shift to the scattering length a and the effective range tq is given by^^ 

k cot 5{k) = -^ + ^rok^ - Prlk^ + Qrlk"^ + ■■■ (2) 

Potentials of different shapes may have identical values of a and tq, but the parameters P 
and Q are shape- dependent. Figure 1 schematically illustrates the zero-energy wave function 
u{r) for an attractive square- well potential, 

V(r) = -K,e(6-r), (3) 

of depth Vo and range 6; Q{x) is the usual unit Heaviside function. 

Figure 1 depicts a situation where the intercept of the the zero-energy wave function (or 
its extrapolation) on the horizontal axis is the scattering length a, with its sign signifying the 
presence or absence of a bound state. For a shallow potential that supports no bound state, 
the scattering length a is negative (Fig. 1(a)). When the depth of the potential is increased 
and fine-tuned as in Fig. 1(b), the system has a zero-energy "bound" state corresponding 
to a ^ ±oo. For this limiting case the dimensionless quantity 7 = MVob"^ /fi^ equals 71^/4. 
The corresponding values of 7 for other two-parameter attractive potentials are listed in 
Ref. Il7|. A further increase in the depth results in supporting only one bound state, and the 



scattering length a is shown to be positive (Fig. 1(c)). 



4 




(c) 



a>0 







r — 






< 


a 



r 



FIG. 1: An attractive spherical square well potential. The thin solid line indicates the zero-energy 
wave function u{r), and the dashed line the extrapolated asymptote, (a) Shallow well with negative 
scattering length, (b) Depth tuned to Vq = 7r^7i^/4M6^ with a ^ oo corresponding to the zero 
energy bound sate, (c) A deeper well supporting only one bound state and positive scattering 
length. Note that u{r) ~ (r — a) for r > b. 

The situation is more complicated when the depth of the potential becomes still stronger 
to support two or more bound states. For the case of the square-well example, it may be 
shown that the s-wave scattering length is given by^ 

/ tan ^/7\ , , 

.^.(l-^), (4, 

Therefore = 3tt/2 signals the appearance of the next zero-energy bound state, with the 
scattering length changing sign. A further increase results in two bound states. Note from 
Eq. (jl]) that the scattering length a depends on the range b of the potential and its strength 
7. In what follows we shall be interested in taking the zero-range limit of the interaction, 
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while keeping the strength parameter 7 fixed. 

The effective range tq is proportional to the range b of the potential, and goes to zero as 
6-^0. For this case Eq. ([2]) reduces to 

k cot 6 = — . (5) 
a 

For positive a we may extrapolate Eq. ([5]) to a bound state. Note that the S-matrix for 
the s-state is defined as S{k) = exp{2iS{k)), and the S-matrix has a simple pole for each 
bound state in the upper half of the imaginary fc-axis on the complex fc-plane.— A zero- 
range potential can support only one bound state. At a bound state therefore, k = in, and 
cot 5 = i. The relation ([5]) for a zero-range potential then reduces to k = 1/a. Hence the 
binding energy B = —E, and the bound state wave function ip{r) is given by 

i,(r) = ?^^^tlM, (,>o). (6) 

r 

B = h^/Ma\ (7) 

Equation ([6]) is exact for a zero-range potential. If tq is nonzero, the relations for B and 
the asymptotic (r > b) wave function are approximately true only if the binding B is 
close to zero, which implies that a ^ tq, and the system is large in size compared to the 
range of the potential. As an example where these conditions are only marginally satisfied, 
consider the deuteron, which is the lightest stable atomic nucleus, and is a bound state 
of a neutron and a proton. The n-p scattering length in the s-state {S = 1, T = 0) is 
a = 5.43 fm, and tq = 1.75 fm, where Ifm = 10~^^m. Equation ([7]) for the binding energy 
gives B = 1.41 Mev, compared with the experimental value of 2.23 Mev.— 

Because tq is nonzero but is much smaller than a, this agreement could be improved by 
taking the next term in the expansion in Eq. ([2]). Following the same extrapolation for the 
bound state at A; = in, we obtain the quadratic equation 

1 1 2 

K = - + -roK\ (8) 

Ct Zi 

We select only the root that gives k = when a ^ 00 and ro = and obtain 

K=-{\- v/(l - 2ro/a)) . (9) 
Equation ^ for the deuteron gives B = 2.16 MeV, which is close to the experimental value.— 
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We have assumed that the size of the bound two-body system is much larger than the 
range of the potential responsible for the binding. The details of the potential at short 
distances do not manifest themselves significantly in the wave function, and the binding is 
determined by the low energy scattering parameters. It should therefore be possible to design 
an effective potential that reproduces the shape-independent results we have obtained. Such 
an effective one-parameter pseudopotential is easily derived using the wave function iplr) 
given by Eq. ([6]). From the relation V^(l/r) = —ATT6^{r) we obtain — V^(exp(— r/a)/r) = 
— exp(— r/a)/r + 47rexp(— r/a)(5^(r). If we identify exp{—r/a)/r = ip{r), we may write 

d 1 

- V^^A + A^a5^ir)—ri) = — -i/'. (10) 
or 

We see that the pseudopotential 

Vir) = ^W(r)|-r (11) 

reproduces the earlier results for the wave function and the binding energy of any zero-range 
potential with a given scattering length. In the many-body problem of an ultra-cold dilute 
gas, where k ^ 0, and the average distance between the particles is much larger than the 
range b of the potential, the same regularized pseudopotential in Eq. ffTTl) is widely used.— 
Note that d{rip)/dr in the limit r — is ^/^(O) for a wave function that is well-behaved at 
the origin. In addition, it also gives a finite result for the irregular il)[r) that goes as r~^, 
which is more relevant for our problem. 



B. Two particles in a harmonic potential 

First consider the case of noninteracting particles. Each particle of mass M moves in a 
harmonic potential l/2Mu;^r^. If we make the usual transformations to center-of-mass and 
relative co-ordinates, r = (ri — r2), and R = (ri + r2)/2, and their corresponding canonical 
momenta p and P, we obtain 

where Mem = 2M and fi = M/2. We concentrate on the relative motion and divide the 
relative distance r by ^/h/ fiuj and the energy E by hu. In terms of the dimensionless scaled 
variables x = r/\/2l, where / = ^J%/ [Mu], and i] = 2E/[%lo), the Schrodinger equation for 
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the relative wave function u{r) = r'ip{r) in the s-state is given by 

+ X u = rju. (13) 

The regular ground state solution is uo{x) = xexp(— corresponding to ipi^x) = 
exp(— a;^/2), with r] = 3, or E = jhu. It is easily checked that Eq. f|T3|) is also satis- 
fied for X > by u{x) = exp(— corresponding to ipi^) = exp(— x^/2)/x with 77 = 1 or 
E = ^huj. Although normalizable and lower in energy, it is excluded as a valid solution in 
the noninteracting problem because at r = 0, the kinetic operator acting on 1/r would 
yield a delta function. 

The situation changes for the zero-range interaction that we have been considering. At 
r = 0^ we may ignore the harmonic confinement, and the wave function is given by Eq. ([6]). 
Even for nonzero b the harmonic potential may be ignored if fiuj^b^ <^ huj, which yields 
the condition h/l <^ 1 for the validity of the spectrum derived in the following. For the 
scattering length a — > 00, the wave function is 1/r, and smoothly joins with the irregular 
wave function above for rj = 1 This analytic behavior suggests that the valid (normalized) 
ground state wave function in the limit of a — 00 for r > is the irregular solution. 



(14) 



with E = ^hu. The next excited i = state has one node, and by demanding that 
it be orthogonal to the the ground state given by Eq. (IT^ . may be deduced to be {x"^ — 
1/2) exp(— a:;^/2)/x, with energy E = ^Tioj. There is an infinite tower of such radially excited 
states with £ = 0. This result is confirmed by our numerical calculations with the two-body 
potential 

V{r) = -Vosedi^ir/b). (15) 

This potential has a bound state at = in free space when Vq = 8h'^ /Mr^, and b = ro/2. 
With these parameters, the scattering length a is infinite, and the effective range is tq. 
Figure 2 displays the first five eigenvalues as a function of ro//. Note that as ro — 0, in 
units of hu approaches 1/2, 5/2, 9/2, . . .. 

The tower of states described above, for a 00, may be generalized for any given a, 
following Ref . [8|. From our previous considerations it is reasonable to assume that the general 
solution for x > is of the form u{x) = exp(— x^/2)/(x). It is convenient to let y = x"^, and 
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FIG. 2: The lowest five i = energy levels for the potential V{r) = -(8/i /MrQ)sech"' (^T^j, as 
ro is decreased toward zero. In this and all subsequent figures, lengths and energies are scaled as 
described in the text. 

w{y) = f{x). After some algebra, Eq. (fT3!) is transformed to 



y-n^ + (o - y)-r + i^—)^ = 0- 



(16) 



Equation ( fT6l) is the equation satisfied by the confluent hypergeometric function, and the 
solution is given by^^ w{y) = ciM((l - r])/A,l/2,y) + C2y^/'^M{{3 - 7])/4:,3/2,y), where ci 
and C2 are constants. We transform back to the x co-ordinate and obtain 



u 



(x) = [ciM((l - r])/A, 1/2, x^) + C2xM{{3 - r/)/4, 3/2, x^)] exp(-xV2). (17) 



The ratio Ci/c2 may be determined by noting that for large y, M{p,q,y) 
T{q) /T(p)y^P^''^ exp^y). The large x behavior of the solution u{x) given in Eq. (ITTj) is 



ulx) ~ 



Cl 



r(l/2) 



+ C2 



r(3/2) 



r((l-ry)/4) ^^r((3-r/)/4)J 



a;-(l+'7)/2gxp(3,2/2). 



(18) 
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For u{x) not to diverge for large x, the quantity in the square brackets must vanish, yielding 
the ratio to be 

cx_ i r((l-r^)/4) 

C2 2r((3-r^)/4)- ^ ^ 

We now relate this ratio to the scattering length a to determine the eigenvalue t]. We 
already noted that as r ^ 0^, the oscillator potential goes to zero, and u{r) for positive 
energy is a scattering solution outside the zero-range potential, given by Eq. ([1]). For small 
r we have to within an overall constant, 

ti(r) ~ (r + i tan 5). (20) 

In contrast, because the confluent hypergeometric functions go to unity for x ^ 0^, we find 
from Eq. f|T7|) . u{x) ~ {ci + C2x). If we factor C2 out of this expression and rewrite u in terms 
of the variable r, we find to within an overall constant, 

u{r) ~ ( r + v^/-). (21) 

V C2 / 

We compare Eqs. (!20l) and (!2T1) and obtain 

-tan5 = V2/ — . (22) 

We substitute tan 5 = —ak in Eq. (!22l) and the ratio Ci/c2 from Eq. (fT9l) and obtain the 
desired relation 

a ^ 1 r((l-r^)/4) 

/ v^r((3-r^)/4)- ^ ^ 

The Gamma function diverges when its argument is zero or a negative integer. Hence, for 
a — > oo, rj = 1,5,9, .. ., confirming our earlier results for both the energy spectrum and the 
corresponding wave functions u{x). 

Note that Eq. (|23l) has been obtained with no mention of the shape of the potential 
and is therefore applicable for any short-range two-body potential. This general analysis 
is in contrast with the derivation in Ref. in which explicit use of the properties of the 
pseudopotential (ITT]) was required. 

Figure 3 displays the variation of the eigenvalues E (in units of hu) as a function of a/l. 
For comparison we also give the values obtained numerically for a square-well potential with 
b/l = 0.01. On the scale of the plots, it is difficult to distinguish between the results. Note 
also that the ground state energy diverges toward —oo as a ^ 0. This behavior is the same 
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FIG. 3: The energy spectrum of the two-particle system in a harmonic trap. The sohd hnes are for 
the attractive square-well potential ([3]) with b/l = 0.01, and the dashed lines for the pseudopotential 
(jlip . This theoretical spectrum with varying scattering length is in excellent agreement with the 
experimentally measured binding energies of the dimers formed from ^'^K atoms.— 

as that of the divergence of the single bound state for the regularized delta function 
potential given by Eq. (I6l). Note that a zero scattering length does not necessarily mean 
that the interaction is zero. If the scattering length approaches from the negative side 
(see Fig. 1(a)), the interaction becomes weaker and eventually vanishes. If it approaches 
zero from the positive side (see Fig. 1(c)), the interaction becomes stronger and stronger, 
squeezing the wave function more and more inside the potential. The flow of the eigenvalue 
curves as a function of a are to be understood with these observations in mind. Starting 
from the left in Fig. 3 with large negative a, the interaction becomes weaker a increases, 
eventually reaching zero when a ^ 0^. At this value we recover the unperturbed eigenvalues 
of the harmonic oscillator. To explain the behavior of the curves in Fig. 3 for a > 0, start 
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from the right-hand end and decrease a. We have aheady commented on the diverging 
ground state. The next one, starting at 5/2 at the right end, continually becomes less as a 
decreases, because the interaction becomes stronger, eventually reaching the lower value 3/2 
at a = 0"*". The previous arguments should help to demystify the counterintuitive features 
of the pseudopotential discussed in Ref. ln. 



III. CONCLUSIONS 



We have discussed the two-body problem in the context of ultra-cold atoms. We have 
shown that the energy spectrum is a universal property of the interacting two-body problem 
provided that the range of the interaction is much smaller than the oscillator length. We have 
clarified the relation between the irregular solution and the regularized 5-pseudopotential. 
We invite the interested reader to apply the methods used here to the interesting one- and 
two-dimensional analogues of this problem. The energy spectra for these cases are also given 
0. 



m Ref. M. The two-dimensional case is elaborated in Ref. 



221 . We suggest Ref. 



23 



and 



24 for 



background reading on one- and two-dimensional scattering. 
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